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o : 

^ ■ Abstract 

> : 

O ■ Second order SUSY transformations between real and complex potentials for three im- 

. portant from physical point of view Sturm-Liouville problems, namely, problems with the 

^3 I Dirichlet boundary conditions for a finite interval, for a half axis and for the whole real line 

• are analyzed. For every problem conditions on transformation functions are formulated when 

. transformations are irreducible. 

> : 

^ 1 Introduction 

(N : 
o . 

^ I Non-Hermitian operators started to attract attention of physicists soon after the main foundation 

■ of quantum theory was built (see e.g. [HE])- A recent numerical observation [3] that some non- 

p ■ Hermitian one- dimensional Hamiltonians may have purely real spectrum re-initiated an attempt 

Jj . to generalize quantum mechanics by accepting non-Hermitian operators for describing physical 

^ I observables [1] (so called 'complex quantum mechanics', for recent developments see e.g. [316]). 

^ I In this respect it is worthy of special mention the paper [7] where the authors establish a general 
criterion for a set of non-Hermitian operators (so called quasi-Hermitian) to constitute a consistent 

^ ■ quantum mechanical system with a normal quantum mechanical interpretation. 

^ I Another class of non-Hermitian operators, called pseudo- Hermitian, was introduced by Dirac 

' and Pauli and later used by Wick and Lee [T] to overcome some difficulties related with using 
Hilbert spaces with an indefinite metric. Recently due to Solombrino and Scolarici ^ this concept 
found a further generalization as weak pseudo-Hermiticity. A recent observation [9] that the real 
character of the spectrum of a pseudo-Hermitian Hamiltonian h is related to the existence of a 
pseudo-canonical transformation, which makes h similar to a Hermitian operator, permits us to 
suppose that there exists a certain overlap between these two classes of non-Hermitian operators, 
appearing to be the most appropriate candidates for describing physical observables. 

From the general point of view if some non-Hermitian operators with a purely real spectrum 
are similar to Hermitian ones their incorporation into quantum mechanics cannot be considered as 
a more general approach with respect to the conventional quantum mechanical description. Their 
use may give (or may not give, for a recent discussion see [10]) only calculational advantages [3 [7]. 
From the first glance this observation leads to the negative answer to the question whether or not 
complex quantum mechanics is an extension of the conventional quantum mechanics. But if we 
take into account the fact that between non-Hermitian operators with a purely real spectrum there 
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exist such operators which never can be similar to Hermitian ones the above question seems to be 
still open. In particular, a non-diagonalizable operator has an incomplete system of eigenfunctions 
and therefore it can never be similar to a Hermitian operator, the eigenfunctions of which form 
a complete basis in corresponding Hilbert space. Other operators which should be studied in 
this respect are the ones having spectral singularities in the continuous part of the spectrum, the 
feature which never appears in the Hermitian case [TT] . 

In [12] an overlap between PT-symmetric quantum mechanics (by this term some authors mean 
a complex extension of quantum mechanics) and supersymmetric quantum mechanics (SUSY QM, 
for recent reviews see [T3]) was noticed. Furthermore, as it was shown SUSY QM may be useful 
to transform non-diagonalizable Hamiltonians into diagonalizable forms [13] and remove spectral 
singularities from their continuous spectrum [15]. These results permit us to hope that SUSY QM 
may be a very useful tool in complex quantum mechanics. 

In this letter we are using the term SUSY transformations in its narrow sense as differential 
transformations between two (exactly solvable) Hamiltonians having almost the same (up to a 
finite number of levels) spectra and do not discuss underlying algebraic constructions. 

Supersymmetric transformations involving first order intertwining operators between one- 
dimensional Hamiltonians 

ho^i = ~dl + Vq^i{x) , = d/dx , a; G (a, h) (1) 

with the possibility for the potentials Vo,i(x) to be complex- valued functions were studied in 
[T6| [T7j. A succession of two SUSY transformations with equal factorization constants (confluent 
transformations, see e.g. [IB]) was used in [IS] to obtain bound states embedded into continuum 
of scattering states of a complex potential. On the other hand it is clear that if the intermediate 
Hamiltonian /i of a chain of two first order transformations ho ^ h ^ hi has a physical meaning 
there are no special needs to study the second order transformation leading from ho to hi directly 
(2-SUSY transformation); all properties of the Hamiltonian hi can be understood at the level of the 
first order transformation h ^ hi (1-SUSY transformation). This is not the case if h = —d'^ + V{x) 
is not a well defined Hamiltonian acting in the same Hilbert space as ho and hi. In this respect for 
the case when both Vo{x) and Vi{x) are real in [20] the notion of reducible and irreducible SUSY 
transformations was introduced. The chain is called reducible if V{x) is real [20] and irreducible 
otherwise. Evidently as far as complex potentials are concerned chains irreducible in this sense 
become reducible [16]. Later [21] (for a recent discussion see [22]) another type of irreducible 
transformations was described. They appear when the potential V{x) has singularities inside the 
interval {a,b) where the potentials Vo{x) and Vi{x) are regular. Recently the third possibility 
for irreducible chains was noticed [23]. It appears if the intermediate potential V{x) is regular 
inside (a, b) but the spectrum of h is completely different of the spectrum of ho and no SUSY 
transformations between their eigenfunctions exist whereas ho and hi are (almost) isospectral and 
their eigenfunctions are connected with the help of a second order SUSY transformation. For 
these two kinds of irreducible chains spectral properties of hi cannot be derived from spectral 
properties of h even in the case of complex potentials and one needs to analyze second order 
SUSY transformations between ho and hi without involving the intermediate potential V{x). In 
this Letter we formulate conditions for 2-SUSY transformations between a real potential Vq and 
a complex potential Vi with a purely real spectrum to be irreducible. 
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2 First and second order SUSY transformations 



In this section we briefly review the main properties of first and second order SUSY transformations 
(for details see e.g. |i3l UHl El] ) we need below. 

We consider second order ordinary differential equations (Schrodinger equations) 

ih-E)fEix)=0 xe{a,b) (2) 

with (Hamiltonian) h = ho,h,hi and = iPe{x),4'e{x),ipe{x) respectively; G C is a 

parameter and a, b may be both finite and infinite. 

We say that the Hamiltonians ho and h are 1-SUSY partners if there exists a first order 
differential operator Li intertwining ho and h, Liho = hLi. Similarly, h and hi are 1-SUSY 
partners if there exists a first order differential operator Li such that Lih = hiLi. Evidently, the 
second order differential operator L = LiLi intertwines ho and hi 

Lho = hiL (3) 

and hi is called 2-SUSY partner for ho. Once the existence of L is established solutions of equation 
([2]) with h = hi can be found by applying L to solutions of the same equation with h = ho, 
ifE = LipE, i'E ^ kerL. Evidently, similar property takes place for solutions ipE of the Schrodinger 
equation with the intermediate Hamiltonian h\ they are expressed in terms of solutions of the 
initial equation ipE and hoUi{x) = aiUi{x), ai G C [ui is called transformation function 

and ai is factorization constant) 

'iPe{x) = LiipEix) = -iPe{^) + w{x)'iI)e{x) w{x) = Ui{x)/ui{x) E ai (4) 



Ui{x) ' 

The potential V{x) is expressed in terms of the function w{x) as follows: 

V{x) = Vo{x) - 2w'{x) . (5) 



For the next transformation step h ^ hi the same formulas (j4]) and ([5]) with the evident 
modifications should be used with the only difference that now one distinguishes the confluent 
case, when the factorization constant a2 at the second step of transformation coincides with that 
of the first step, a2 = oii, from the usual (non-confiuent) case when these constants are different 

Using the second order transformation operator L = LiLi one can avoid the intermediate step 
and go from ho to hi directly 

Vi = Vo-2[\ogWiui,U2)]" (6) 

(fE = Li)E = W{ui,U2,i'E)/W{ui,U2) . (7) 

Here and in the following the symbol W with arguments being functions denotes Wronskians, 
ui, U2 and ipE are solutions to equation ([2]) with h = ho corresponding to the eigenvalues ai, 0.2 
(factorization constants), /io'"i,2 = «i, 2^*1,25 and E respectively. Expressions ([6]) and ([7]) are known 
as particular cases of Crum-Krein formulas |25j . 
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Formula ([7j) defines the operator L for any sufficiently smooth function ipE but if is a 
solution to equation ([2]) other forms of this equation are useful 

<^E = - "2) V'B + ("1 - "2)7777 tU\ (8) 

= [E - aijipE + (ai - "2) 7777 r M2 • (9) 

W[Ui,U2> 

Here the use of equation ([2]) has been made to express the second derivatives of the functions 
Ui, U2 and V'e ill terms of the functions themselves. Operator L as given in ([8}l9]) maps any two- 
dimensional space of solutions of equation ([2]) with h = ho at E ^ ai, a2 onto corresponding space 
of solutions of the same equation with h = hi. The two-dimensional space span(ui,M2) is the 
kernel of L, Lui^2 = 0. Despite that with the help of L one can find solutions of the transformed 
equation corresponding to E = ai,a2- For this purpose one has to act with L on functions 
Vi,2 7^ Ml, 2, hQVi^2 = "1,2^^1,2- Using the fact that 1^(^1,2,^1,2) = const and putting ipE = "^1,2, 
E = ai,2 in and ([9]) one readily gets 

where we have omitted an inessential constant factor. It is worth to note that the use of these 
functions for the next step of the second order transformation gives back the initial Hamiltonian 
ho and, hence, the procedure is completely reversible. Our last comment here is that as it follows 
from ([6]) to obtain nonsingular for x G (a, b) potential differences it is necessary that W{ui, U2) 7^ 
Vx G (a, h) which will be supposed to be the case. 

For equal factorization constants (ai = 02 = ", the confluent case) the function W{ui^U2) in 
(E]) should be replaced by 

W,{x) = c+ I u\y)dy. (11) 



Constants c and Xq should be chosen such that Wc{x) 7^ Vx G (a, h). Solution of the Schrodinger 
equation with the Hamiltonian hi are given by 

i^E{x) = L^E = {a-E) tljEix) + ^^^^^^ ^i^) ■ (12) 

This formula gives a solution for ii^ = a also provided ipa is linearly independent with u 

^.{x) = (13) 
Wc{x) 

where once again we have omitted an inessential constant factor. 

The properties described above take place irrespective of any boundary value problem related 
to the differential equation ([2]). Here we shall consider boundary value problems of three kinds 
which are the most interesting from physical viewpoint: 

(i) Regular Sturm-Liouville problem; the potential Vq is bounded and continuous in [a, b] which is 
a finite interval. We will consider only Dirichlet boundary conditions, iPe{o) = 0, ipEib) = 
imposed on smooth (infinitely differentiable) functions from L'^{a,b) which form an initial 
domain of definition of ho which has a purely discrete spectrum. 
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(ii) Scattering potentials on a semiaxis, i.e. Vo{x) is continuous and bounded from below for 
X G [0, oo) and such that 



Here we impose on solutions to equation ([2]) the Dirichlet boundary condition at the origin 
only, iPe{0) = 0, which together with the condition of square integrability over the interval 
[0, oo) selects the bound states. The scattering states have an oscillating asymptotical be- 
havior. The initial domain of definition of ho consists of infinitely differentiable functions 
vanishing for sufficiently large x and at the origin. The operator ho has a finite number of 
discrete levels and the continuous spectrum filling the positive semiaxis. 

(ill) Confining and scattering potentials on the whole real line, (a, b) = M. For confining potentials 
Vq is locally bounded and Vo(x) — > cxd when — > oo. Scattering potentials are selected by 
the condition 



where Vq{x) is continuous and semi-bounded from below. Operator is initially defined on 
the set of infinitely differentiable functions from L^(]R) vanishing for sufficiently large \x\. 
For confining potentials the spectrum is purely discrete. The scattering potentials have a 
finite number of discrete levels and a two-fold degenerate continuous spectrum filling the 
positive semiaxis. 

In all cases the operator is essentially self-adjoint and has a complete set of eigenfunctions 
(in the sense of generalized functions for continuous spectrum eigenfunctions) in corresponding 
Hilbert space (see e.g. [22]). 

3 Complex SUSY partners of real potentials 
3.1 General remarks 

As we shall see both the analysis and the results strongly depend on the character of the initial 
Sturm- Liouville problem. Nevertheless, there is a property common to all eigenvalue problems 
which is essential for our analysis. We formulate it as the following 

Proposition 1. For a real potential Vo{x) defining a self-adjoint operator Hq in the space L'^{a, b) 
any solution V'a(x) to the Schrddinger equation with Im (a) ^ can vanish at no more then in 
one point of the interval [a,b]. For nonfinite values of a,b the statement should be understood in 
the sense of limit, i.e. for instance for b = oo this is limx^oo ipaix) = 0. //Im(Q!) = and exists 
Xq e [a,b] such that ipaixo) = then up to an inessential constant factor ipa{x) is realWx G [a,b]. 

Proof. The first part of the statement follows from the property of a self-adjoint operator to 
have only real eigenvalues. Indeed, if equations V'a(xo) = ipaixi) = took place for Im(a;) ^ 
and xo,xi e [a,b], then the operator ha, defined by the same differential expression ho and the 
zero boundary conditions at the ends of the interval [a;o,a;i], being self adjoint would have the 
complex eigenvalue E = a which is impossible. 




(14) 




(15) 
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The second part follows from the property that for a real a the basis functions, il)ia{x) and 
ip2a{x) in the 2-dimensional space of solutions of equation ([2]) can always be chosen real so that 
any complex- valued solution ^ipaix) is a linear combination ^ipaix) = ciipia{x) + C2ip2a{,x). For a 
finite value of Xq from the equation ipa^Xo) = one of the constants, say C2 (if i'2a{xo) 7^ 0) can be 
found. Evidently, it is proportional to Ci with a real proportionality coefficient. If ipai^oo) = 
the function ipa{x) is real- valued for a real a up to an inessential constant factor. The statement 
follows from a contradiction which appears if one supposes the opposite statement to be true. 
Indeed, if ipa{x) is a complex- valued function then ip^i^) (asterisk means complex conjugation) is 
linearly independent with ipa{x) and ?/'*(±cxd) = which is impossible. □ 

We impose on solutions to equation ([2]) the same zero (Dirichlet) boundary conditions after 
2-SUSY transformation. Thus we have two boundary value problems, initial and transformed, 
which we will denote (I) and (II) respectively. 

For the usual (non-confluent) case, 2-SUSY transformation is reducible if 1-SUSY transforma- 
tion with either u = ui ot u = U2 is 'good'. This is due to the fact that the chain can start from 
either u = Ui or u = U2- This contrasts with the confluent case since the chain starts now always 
from 1-SUSY transformation based on the transformation function u. Therefore any 2-SUSY 
transformation is irreducible if this transformation is 'bad'. 

Below we analyze conditions for the transformation functions Mi 2 and factorization constants 
ai^2 giving (according to (jH])) for a given real Vq{x) a complex potential function Vi{x) such that 
the operator hi defined in the corresponding Hilbert space according to the cases (i)-(iii) has a 
real spectrum coinciding with the spectrum of ho with the possible exception of one or two levels 
and the transformation is irreducible. We illustrate every possibility with the simplest example 
considering boundary value problems for Vo(x) = 0. 

3.2 Regular Sturm-Liouville problem 

The problem (I) is regular but the problem (II) may become singular if the potential Vi{x) is 
unbounded in one of the bounds of the (finite) interval [a, h] or in the both. 

To distinguish irreducible second order transformations from reducible ones we have to start 
the analysis from first order transformations. From it follows that if (a) =0 then iPe{0') = 
if and only if Ui{a) = 0. Hence to keep the zero boundary conditions we have to choose the function 
Ui{x) vanishing both at a; = a and at x = b. This means that it is an eigenfunction of the problem 
(I). Since any eigenfunction except the ground state function has zeros inside the interval {a,b) 
we conclude that in this case there exists the only admissible first order SUSY transformation. 
It corresponds to Ui{x) = ipo{x) (the ground state function) after which the ground state level is 
deleted. This means that any 2-SUSY transformation which does not involve the ground state 
function of the problem (I) is irreducible. 

It is clear from (^M> that if both ui (or equivalently U2) and ipE, E ^ ai, 0:2, satisfy the zero 
boundary conditions then (fE satisfies the zero boundary conditions also. Hence, to keep the zero 
boundary conditions after 2-SUSY transformation we have the following possible choices for ui 
and U2: 

(a) ui{a) = ui{h) = U2{a) = U2{h) = 0; 

(b) ui{a) = ui{h) = 0, U2{a) = (or ^2(6) = 0), U2{h) ^ (or ^2(0) ^ 0); 
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(c) ui{a) = ui{h) = 0, U2{a) ^ 0, U2{b) ^ 0; 

(d) ui{a) = U2{h) = (or ui{b) = U2{a) = 0), ui{b) ^ 0, ^2(0) ^ (or ui{a) ^ 0, U2{b) ^ 0). 



In case (a) both ui and U2 are eigenfunctions of the problem (I) and there is no way to get a 
complex potential difference. 

In case (b) ui is an eigenfunction of ho, ui{x) = ipk{x) and ai = Ek, k = 0,1, . . .. Hence, the 
level Ek is not present in the spectrum of hi. The parameter 02 can take any complex value except 
for a2 = El, I = 1,2 I k since in this case U2 = ipi and we are back in the conditions of the 
case (a). According to ffTOj) the function ipa^ satisfies the zero boundary conditions and, hence, the 
point E = a2 belongs to the discrete spectrum of hi. Using Proposition 1 we conclude that in this 
case a complex potential Vi is possible only for a complex value of 02. The Hamiltonian hi has 
thus the complex discrete level E = 02. So, there are no ways to obtain a complex potential with 
a real spectrum in this case. The 2-SUSY transformation is reducible if Ui{x) = ipolx) (ground 
state function) and irreducible otherwise. 

In case (c) ui is still an eigenfunction of ho, ui{x) = ipk{,x), ai = Ek, k = 0,1, . . . and the level 
Ek is not present in the spectrum of hi but there are no restrictions on U2. Yet, the level E = a2 
belongs to the spectrum of hi and if we want for the Hamiltonian hi to have a real spectrum we 
have to choose a2 real. In this case complex potential differences can arise from formula ([6]) only 
if U2 is a complex linear combination of two real linearly independent solutions of equation ([2]) 
with h = ho. As it is shown in [T3] if C(2 = Ei the Hamiltonian hi becomes non-diagonalizable. 
The second order transformation is irreducible provided A; > 0. Another interesting feature we 
would like to mention is the possibility to get PT-symmetric potentials by appropriate choice of 
the function U2{x) and if Vo(x) has this property. 

Example 1. For x G [— vr, tt] take Ui = sin(noa;), (cti = ng). uq = 1,2, . . ., U2 = cos(ax + b) 
(02 = a^), a G M, a 7^ uq, Im (6) 7^ 0. The potential Vi is given by 

2 2\ nl\cos{2ax + 2b) + 1] + a'^\cos{2nox) - 1] 

^1 = K - « ) 7 7 ^ -f -TT— —f T— 7 -TTT^ • (16) 

[rio cos(noa;j cos(ax + 0) + a sm(noxj sm(ax + b)\'^ 

If a 7^ ±n/2, n = 1, 2, . . . the spectrum of the Hamiltonian hi with potential fll6p consists of all 
levels of ho = —d"^ which are -En-i = except for E = and an additional level Ea = a?. If 
Re (6) = this potential is explicitly PT-symmetric. 

One can find other examples of potentials one can get under these conditions in [T3|. 

Consider finally the case (d). Using Proposition 1 we conclude that if both ai and 02 are real 
there is no way to obtain a complex potential difference. So, to be able to produce a complex 
potential Vi we have to choose at least one of as (say ai) complex. In this case Ui is nodeless and 
the intermediate Hamiltonian h is well-defined in L^(a, b) but its spectrum is completely different 
from the spectrum of /iq since corresponding 1-SUSY transformation breaks the zero boundary 
condition aX x = b. We thus can construct an irreducible SUSY model of a new type. Here also 
one can get "PT-symmetric potentials if Vo{x) is PT-symmetric and a2 = «*■ This is readily seen 
for a symmetric interval b = —a. The last comment here is that according to ffTOj) neither ipa^^ nor 
(Pa2 satisfy the zero boundary conditions. Therefore the Hamiltonian hi is strictly isospectral to 
ho and, hence, its spectrum is purely real. 
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Example 2. For x G [— vr, vr] taking ui = sin(ai(x + vr)) (ai = al) and U2 = sin(a2(a; — vr)), 
{a2 = al), ai 7^ 02, Im (af ) 7^ 0, Im (02) 7^ one gets the following potential: 




a^l — cos(2ai(x + vr))] — af[l — cos(2a2(x — tt))] 



(17) 



[ai cos(ai(x + n)) sin(a2(x — n)) — 02 cos(a2(x — vr)) sin(ai(x + vr))]^ 



which is PT-symmetric provided 02 = a^. The spectrum of the Hamiltonian hi coincides with 
the spectrum of ho, En-i = n'^/A, n = 1,2,.... Although 1-SUSY transformations both with 
u = Ui and u = U2 produce potentials regular in the interval (— vr, vr) they do not preserve the zero 
boundary conditions at both limiting points of (a, h) and therefore the 2-SUSY transformation is 
irreducible. 

For the confluent case as it follows from ( fT3l) the function u{x) should vanish at either bound 
of the interval [a, h] and therefore it is one of the eigenfunctions oi h^, u = ipk and a = E^, 
k = 0, 1, . . .. Therefore to obtain a complex potential Vi ([6]) one has to choose the constant c 
in f|TT]) complex. This transformation is irreducible provided u ^ ipQ. It keeps the spectrum 
unchanged since the function ( fT3l) satisfies the zero boundary conditions. In some cases (i.e. 
Vq{—x) = Vo(x), b = —a, xo = and Re(c) = 0) the potential Vi{x) is PT-symmetric. 

Example 3. For x G [— vr, vr] taking Xq = we get the potentials 



Here the upper sign corresponds to u = cos(nox/2), no = 3, 5, 7, . . . and the lower sign corresponds 
to u = sin(nox/2), uq = 4, 6,8, . . ., (a = nl/A), Im (c) 7^ 0). For Re(c) = these potentials are 
PT-symmetric. The Hamiltonian hi is isospectral with ho = —d^. 

3.3 Scattering potentials on a semiaxis 

As it was already mentioned at the beginning of Section 3.2, 1-SUSY transformation keeps un- 
changed the zero boundary condition (at the origin in the current case) if Mi(0) = 0, Ui{x) 7^ 
Vx G (0, 00). The zero boundary condition at the infinity for a transformed function is satisfied for 
any 1-SUSY transformation provided the transformation operator acts on a function vanishing at 
the infinity. Therefore any 2-SUSY transformation involving a transformation function vanishing 
at the origin and nodeless in the positive semiaxis is reducible. 

Consider spectral problem (II). To keep the boundary condition at the origin after 2-SUSY 
transformation according to ([H]) one has to impose the same condition on one of the transformation 
functions, say Ui{x), i.e. mi(0) = 0. For the second order transformation to be irreducible 
one has to take care of presence of a positive node in ui{x). According to Proposition 1 if 
«! is complex the function ui{x) is nodeless in (0, 00) and 2-SUSY transformation is reducible. 
Therefore to construct an irreducible 2-SUSY transformation we have to choose only real values 
for ai. According to the second part of the same proposition the function ui{x) can be chosen 
real without any loss of generality. So, we choose ai real, Ui{x) real- valued and for ^2 7^ «i we 
enumerate the following possible choices for U2'- 



Vi = TUo 



2 no(2c -|- x) sm{nox) + 2 cos(r;,ox) ± 2 



(18) 



[sm{nQx) ± no (2c + x)]'^ 



(a) U2{0) = 0; 

(b) ^2(0) 7^ and M2(oo) = 0; 
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(c) ^2(0) 7^ and U2{oo) = 00; 

(d) ^2(0) 7^ and U2{x) has an oscillating asymptotics at the infinity. 

For all cases (a)-(d) if ai is a point of the discrete spectrum of ho, the function ui is a 
eigenfunction of the problem (I) (it is a bound state) and the point E = ai does not belong to 
the discrete spectrum of the problem (II). 

In case (a) by the same reason as it was explained above an irreducible 2-SUSY transformation 
is possible if ^2 is real but according to Proposition 1 it may produce only a real potential 
Vi- Nevertheless, one can get interesting complex potentials by reducible transformations with a 
complex 0^2. If ai does not belong to the discrete spectrum of Hq (i.e. |ui(cx3)| = 00) the 2-SUSY 
transformation is isospectral. 

Example 4- Choose ui = sm{kox) (ai = /cq > 0) and U2 = sinh(ax) (0:2 = G C, Ima2 7^ 0). 
Formula ([6]) gives the potential 

^ (^2 ^ fcg[cosh(2ax) - 1] + a^[cos(2fcox) - 1] 

° [ko cos(/coa;) sinh(aa;) — a sin(A;oa;) cosh(aa;)]^ 

Formula (fTOj) gives an oscillating solution of the Schrodinger equation with the potential (fT9l) at 
E = but it is irregular at the origin. A regular at the origin solution at E = k"^ is unbounded 
as X ^ 00. 

Similar to the case (a) in case (b) with a real 02 according to Proposition 1 the potential 
Vi remains real. Therefore to obtain a complex Vi one has to choose 0:2 complex. The function 
ipa2{x) (flOjl is vanishing at the origin but increasing at the infinity. So, the 2-SUSY transformation 
creates a complex potential with a real spectrum coinciding with the spectrum of ho with the 
possible exception of the point E = ai (if it belongs to the spectrum of hi). Here 2-SUSY 
transformation is reducible only if ui{x) is nodeless in (0, 00). According to Proposition 1 the 
function U2 with a complex 0^2 is nodeless and produces a "good" intermediate potential but 
the first order transformation operator based on U2 does not transform eigenf unctions of /iq iiito 
eigenf unctions of h. So, the 2-SUSY transformation is irreducible if ai is real, the function Ui{x) 
has a node in (0, 00) and ^2 is complex. 

Example 5. Once again we choose Ui = sin(/cox) (ai = /cq > 0) but U2 = e"^ (^2 = — a^, 
Re(a) < 0, Im (a^) 7^ 0). The potential Vi reads 

^ 2fco(fco + g^) .^QN 

[koCOs{kQx) — a sin(A;ox)]^ 

In contrast to the case (b) in case (c) function (fa2 (fTOl) is an eigenfunction of hi and E = a2 
is its spectral point. Therefore to get Hamiltonian hi with a real spectrum we have to choose 
real values for a2- Complex potentials can arise in this case if the function U2 is a complex linear 
combination of two real linearly independent solutions to equation at E = a2- So, the 2-SUSY 
transformation creates a new energy level E = a2- It is reducible if Ui{x) is nodeless in (0, cxd) 
and irreducible otherwise. 

Example 6. Choice Ui = sin(A;oa;) {ai = k^) and U2 = cosh(ax + c) (02 = "O^), kQ,a &W and 
Im(c) 7^ results in the potential 

Y^ = (^^2 _^ ^2^ a^[l - cos(2A;oa:)] + + cosh(2ax + 2c)] ^^^^ 
° [ko cos(A;ox) cosh(ax + c) — a sin(A;oa;) sinh(aa; + c)]^ 
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having a discrete level E = —a?. 

In case (d) a2 > 0. Therefore like in the previous case to have a complex Vi U2 should be a 
complex linear combination of two linearly independent solutions to equation ([2]). The function 
-002 belongs to the continuous spectrum of hi which has a purely real spectrum. In some cases 
indicated in [15J the point E = a2 is a spectral singularity of the Hamiltonian hi. 

Example 7. To illustrate this case we take Ui = sm{kox) (ai = /cq) and U2 = sin(A;ix + c) 
(q!2 = kf), ko, ki G M, fco 7^ ki, Im (c) ^ thus getting the potential 

= _ kfjl - cos(2A;ox)] - A;g[l - cos(2A;i + 2c)] 

^ ° [/co cos(/i;ox) sin(/cix + c) — A;i sin(A;oa;) cos(A;iX + c)]2 



For the confluent transformation ([61), (ITTI) . (1131) to keep the zero boundary conditions we have 
to choose u{0) = 0. Therefore by the same reason as it was explained above to get an irreducible 
2-SUSY transformation we have to choose real values for a leading to a real- valued function u{x). 
To obtain a complex potential Vi one has to choose for the constant c from flTTl) a complex value. 
If u{x) decreases at the infinity (i.e. it is an eigenfunction of hi) the function Wc{x) flTTl) is finite at 
the infinity and ipa ([13]) is an eigenfunction of hi. For xq one can choose both Xq = and Xq = oo. 
Since u{x) is square integrable different choices for xq affect only the value of c. If u{x) increases at 
the infinity, one can choose Xq = 0. Therefore the function Wc{x) increases as when x oo 

and the function (1131) is an eigenfunction of hi too. If a > the function u{x) oscillates at the 
infinity as a linear combination ci exp{—^^/ax) + C2exp{^^/ax). The integrand in (ITTl) increases 
as a linear function of x and the potential Vi keeping its oscillating behavior decreases like l/x"^. 
Therefore it does not satisfy condition (IT^ and it is not a scattering potential. This leads to the 
existence of a discrete level embedded into the continuous spectrum since the function (lT3l) is 
vanishing at the origin and square integrable for x G [0, cxo). So, to get a complex potential by 
confluent transformation ([6]), ( ITTl) one has to choose a real, c complex and xq = 0. If a is not a 
discrete spectrum level then it appears as a new energy level for the Hamiltonian hi. If u{x) is 
nodeless in (0, oo) the transformation is reducible and irreducible otherwise. 

Example 8. We choose u = sin(A;oa;) {a = k^ > 0) and replace c — c/2 (Im (c) 7^ 0). The 
potential 

V. - 32k' sm(k.x) Mkox)-ko{x + c)cosikox) 
Vi - 62k, sm(fcox) ^^.^(2fcox) - 2fco(x + c)]^ 

has the discrete level E = k^ embedded into continuum of the scattering states. 



3.4 Scattering and confining potentials on the whole real line 

For a scattering potential the logarithmic derivative of any decreasing or increasing at x — ±00 
solution u{x) of the Schrodinger equation is asymptotically constant. For a confining potential the 
logarithmic derivative of any similar solution u{x) is usually such that the product ip e{x)u' {x) / u{x) 
increases or decreases together with iPe{,x) so that in the last case its asymptotics is square inte- 
grable; the behavior we will assume to take place. For a scattering potential if u{x) has an oscil- 
lating asymptotics and iIje{,x) has an exponentially decreasing one the product ip e{,x)u' {x) / u{x) 
is exponentially decreasing also so that LijjE belongs to the discrete spectrum of h provided iIje 
belongs to the discrete spectrum of h,. This means that if Ui{x) is nodeless in M it produces a 
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"good" intermediate Hamiltonian h. Thus, if a is real according to Proposition 1 any essentially 
complex-valued solution of equation ([2]) (i.e. a solution u{x) such that u'{x)/u{x) is a complex- 
valued function) is suitable for getting a "good" complex first order potential difference. Therefore 
if both «! and a2 are real, any 2-SUSY transformation which may produce a complex potential 
Vi is reducible. Let at least one of the factorization constants, say ai be complex (another con- 
stant, a2, may be both real and complex). If ui{x) vanishes at one of the infinities then according 
to Proposition 1 it does not have real nodes and 2-SUSY transformation is reducible also. If 
|Mi(a;)| — >• oo when |x| oo then according to ffTOl) the potential Vi has a complex eigenvalue 
E = ai. We conclude, hence, that no irreducible 2-SUSY transformations giving a complex poten- 
tial Vi with a real spectrum exist. Of course this does not mean that such transformations cannot 
create complex potentials with a real spectrum but this means that any such a transformation 
can always be presented as a chain of two 'good' 1-SUSY transformations (with a possibility for 
the intermediate potential to be complex). In particular, if Im (0:1,2) 7^ and the functions ui{x), 
U2{x) vanish at different infinities (i.e. for a scattering potential they are two Jost solutions) 
one can get a complex potential isospectral with Hq (hence, its spectrum is purely real). If for 
Vo{x) = Voi~x) in addition 02 = ^1(2;) is PT-symmetric. The last comment here is that 
irreducible 2-SUSY transformations can produce potentials with complex eigenvalues. 

Example 9. Take ui = sinh(ai(x — Xi)) (oi = —af), U2 = sinh(a2(a; — 0:2)) («2 = "^^2); 
Im(a^) 7^ 0, Im(a2) 7^ 0, ai 7^ 02, Xi,X2 G M. The potential 

Y — a2[l - cosh(2ai(x - Xi))] - a^[l - cosh(2a2(a; - X2)] 

^ ^ [a2 cosh(a2(x — X2)) sinh(ai(x — Xi)) — ai cosh(ai(x — Xi)) sinh(a2(x — 0:2))]^ 

has discrete levels aX E = —a\ and E = —a^. If X2 = —Xi and 02 = al it is explicitly VT- 
symmetric. 

Consider finally the confluent case. If Im (a) 7^ and u{x) increases at both infinities, as 
it follows from ( |T3l) the point E = a belongs to the spectrum of hi. If u{x) decreases at one 
of the infinities then according to Proposition 1 the function u{x) has no real nodes and 2- 
SUSY transformation is reducible. Hence, irreducible 2-SUSY transformations creating complex 
potentials Vi with a real spectrum are possible only with real values of a and Im (c) 7^ 0. In this 
case the function (fT3|) has a decreasing asymptotical behavior both for increasing and decreasing 
u{x) as well as for u{x) having an oscillating asymptotical behavior. Therefore the level E = a 
belongs to the discrete spectrum of hi. Simple potentials one can get in this way correspond to 
the choice of an eigenfunction (e.g. a bound state function) of ho as the transformation function 
u{x) since in many cases it is described in terms of elementary functions. We do not illustrate 
these possibilities by examples and refer an interested reader to existing literature where this is 
done [19]. 



4 Concluding remarks and some perspectives 

In summary, in this Letter a careful analysis aimed to distinguish irreducible transformations be- 
tween all 2-SUSY transformations for three important Sturm-Liouville problems, namely, regular 
problem, problem on a half axis and problem on the whole real line is given. We remind that we 
call irreducible those second order transformations for which either (i) the intermediate Hamilto- 
nian of corresponding chain of two transformations is not well defined in the same Hilbert space as 



11 



the initial and final Hamiltonians or (ii) the intermediate Hamiltonian is well defined but its eigen- 
functions cannot be obtained by acting with the intertwining operator either on eigenfunctions of 
the initial Hamiltonian or on those of the final Hamiltonian. It is shown that for the whole real 
line the only possibility for such a transformation to be irreducible corresponds to the confiuent 
case, i.e. to a chain of transformations with coinciding factorization constants. For problems on 
a half line and on a finite interval there are more possibilities. In particular, transformations of 
type (ii) lead to new irreducible SUSY models. 

Using the property of SUSY transformations to provide us with a general solution of the 
Schrodinger equation at any fixed value of the energy one can observe in examples 4-7 an unusual 
property of SUSY transformations and an intriguing phenomenon concerning spectral properties 
of non-Hermitian operators. In the usual practice of SUSY transformations [131 [HI 121] if a 
transformation function corresponds to a spectral point (for the usual non-confiuent case this 
may be only a point of the discrete spectrum) this point is deleted by the transformation. In 
examples 4-7 we used a transformation function corresponding to a point in the continuous part 
of the spectrum of the initial Hamiltonian but in contrast to conventional SUSY transformations 
(i.e. transformations between Hermitian operators) this point now still belongs to the continuous 
spectrum of the transformed problem. This statement follows from a general spectral theorem 
(see e.g. [27]) according to which the spectrum of a closed operator is a closed set and the 
property that if a point is removed from a closed interval of the real axis it is transformed into 
two (semi-)open subintervals. Of course, the operator hi with the initial domain of definition as 
described in Section 2 is not closed but it is closable and its closure coincides with hf^ (see e.g. 
[28] , theorem 6.3.2). Since the eigenfunctions of the operator coincide with complex conjugate 
eigenfunctions of hi, the operator h^^ has the same system of eigenfunctions as hi. Moreover, the 
analysis of solutions of the Schrodinger equation a.t E = > shows that a solution vanishing 
at the origin is unbounded. This contrasts with the usual quantum mechanical requirement that 
continuous spectrum eigenfunctions should be bounded. We think that instead of imposing the 
above quantum mechanical requirement on continuous spectrum eigenfunctions one should use 
the fact that these functions being ordinary locally integrable functions are in fact generalized 
eigenfunctions of the hamiltonian hi and should be considered as functionals over the domain 
of definition of hi. (In other words one has to involve the notion of the Gelfand triplet into 
analysis [29].) From this point of view the question whether continuous spectrum eigenfunctions 
are bounded or not has no importance. Nevertheless, if necessary one can analyze the growth of 
generalized eigenfunctions for x — > cxd (see theorem 6 in section 55 of [27J) as they are ordinary 
locally integrable functions. 

We think that the results of the present Letter are important in view of the notion of AT-fold 
supersymmetry [30] . As far as this notion is applied to Hamiltonians of type ([1]) with supercharges 
built of differential intertwining operators acting in the same Hilbert space as the components of 
the super-Hamiltonian one can always apply a theorem [21] to factorize an A^th order in deriva- 
tive intertwining operator to a superposition of first order operators thus replacing an A^th order 
transformation by a chain of only first order transformations. If all intermediate Hamiltonians 
of the chain are defined in the same Hilbert space as the initial and final Hamiltonians and at 
any step of transformations the eigenfunctions of two neighbor Hamiltonians are connected by 
corresponding transformation (i.e. intertwining) operator the AA-fold supersymmetry is reducible 
and, actually, all properties of the final Hamiltonian (as well as any intermediate Hamiltonian) can 
be understood at the level of a chain of simpler first order (i.e. usual) supersymmetry transfor- 
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mations. Evidently this is not the case if at least one Hamiltonian of the chain is not well-defined 
(case (i) of irreducible supersymmetry) or at least for one of the Hamiltonians we will not be able 
to get eigenf unctions by applying transformation operator to eigenf unctions of its SUSY part- 
ner (case (ii) of irreducible supersymmetry). The main result of the present Letter consists in 
formulating conditions on transformation functions to produce the simplest irreducible two-fold 
supersymmetry between real and complex potentials. 

Another field of application of our results is the supersymmetric approach to the inverse scatter- 
ing problem (so called supersymmetric inversion [SI])- Since irreducible chains of transformations 
proved to be very efficient in supersymmetric inversion for usual Hermitian operators [32] we hope 
that our results open a way for wider applications of the method of supersymmetric inversion to 
complex potentials previously used for obtaining complex optical potentials only at the level of 
reducible transformations [33] . 
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